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This  report  will  describe  the  Hydrologic  Engineering  Center's  experience  with 
land  use  classification  from  LANDSAT  niul tispectral  imagery.  Land  use  is 
required  for  the  estimation  of  hydrologic  model  parameters.  The  land  use 
classification  procedure  used,  developed  at  the  University  of  California, 
Davis,  for  the  Corps  of  Engineers,  is  an  unsupervised,  noninteractive  approach 
requiring  no  special  image  processing  equipment.  Watershed  land  use  was 
determined  from  LANDSAT  digital  data,  entered  into  a  geographic  data  bank, 
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I  significance  of  accuracy  in  land  use  identi f ication  to  be  assessed. 
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ABSTRACT 


This  report  will  describe  the  Hydrologic  Engineering  Center's 
experience  with  land  use  classification  from  LANDSAT  mul t i spec tra  1  imagery. 
Land  use  is  required  for  the  estimation  of  hydrologic  model  parameters.  The 
land  use  classification  procedure  used,  developed  at  the  University  of 
California,  Davis,  for  the  Corps  of  Engineers,  is  an  unsupervised, 
noninteractive  approach  requiring  no  special  image  processing  equipment. 
Watershed  land  use  was  determined  from  LANDSAT  digital  data,  entered  into  a 
geographic  data  bank,  and  compared  with  a  conventional  land  use 
c 1  ass L f icat ion.  Hydrologic  simulation  model  parameters  were  estimated  from 
land  use  an-J  other  basin  characteristics.  The  generated  discharge  frequency 
curves,  corresponding  to  the  alternative  land  use  classifications,  permitted 
the  hydrologic  significance  of  accuracy  in  land  use  identification  to  be 
assessed. 


( KEY  TERMS :  Land  use ;  LANDSAT ;  hydrologic  simulation;  spatial  data 
management . ) 
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INTRODUCTION 


The  hydrologic  modeling  of  a  watershed,  particularly  urban  or  urban¬ 
izing  basins,  requires  that  the  distribution  of  land  use  be  determined.  The 
amount  and  timing  of  runoff  is  directly  related  to  the  infiltration  capacity 
of  a  land  area  with  the  most  important  distinction  being  between  pervious 
and  impervuous  land  surfaces.  Water  quality  parameters  have  a  similar  df*- 
pendence  on  land  use  data;  rate  of  accumulation  of  a  particular  pollutant 
per  unit  area  is  normally  expressed  as  a  function  of  land  use.  Water  re¬ 
source  planning  studies  are  interested  in  not  only  an  assessment  of  the 
present  state  of  the  water  and  related  resource  system,  but  also  its  pos¬ 
sible  future  configuration.  By  expressing  hydrologic  parameters  as  a 
function  of  current  land  use  it  becomes  possible  to  rationally  predict  the 
impact  future  Land  use  changes  will  have  on  the  quantity  and  quality  of 
future  runoff. 

Manual  methods  for  land  use  ident i f icat i on  Ce.g.,  interpretation  of  low 
altitude  aerial  photography  and  field  surveys)  are  frequentlv  used  in  water¬ 
shed  studies.  The  problem  with  this  approach  is  that  the  resource  re¬ 
quirements,  both  money  and  labor,  for  manual  classification  can  be  exten¬ 
sive.  An  attractive  alternative  is  the  utilization  of  available  remote  sen¬ 
sing  systems  and  computer-ass i sted  classification  techniques.  The  LANDSAT 
satellites  have  been  shown  to  have  the  capability  of  providing  land  use  data 
at  acceptable  levels  of  accuracy  for  hydrologic  modeling  purposes  (Ragan, 
1975;  Jackson,  1977),  LANDSAT  data  is  quicker  and  less  costly  to  obtain  and 
interpret  than  low  altitude  aerial  photography,  provides  repetitive  coverage 
of  the  same  area  at  least  every  18  days,  and  is  available  for  the  entire 
United  States.  Additionally,  lANDSAT's  digital  format  can  be  directly 
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analyzed  by  the  diff(»rent  classification  computer  programs  availahlr,  an<i 
can  be  resampled  for  automatic  inclusion  in  a  geographic  data  bank. 
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This  paper  describes  the  land  use  classification  procedure  developed 
for  the  Corps  of  Engineers  by  the  University  of  California,  Davis  UCD )  and 
presents  results  from  applying  this  procedure  to  land  use  classifications  of 
a  Texas  watershed.  An  extensive  set  of  computer  programs  developed  by  thr^ 
Hydrologic  Engineering  Center  CHEC)  for  the  automated  analysis  of  hydro¬ 
logic,  economic,  and  environmental  aspects  of  comprehensive  water  resource 
planning  studies  based  on  spatial  data  were  used  to  compare  the  LANDSAT 
classification  to  conventionally  determined  land  use.  The  calculated  dis¬ 
charge  frequency  curves,  corresponding  to  the  LANDSAT  and  conventional  land 
use  classifications,  are  presented  as  a  measure  of  the  hydrologic  signi¬ 
ficance  of  differences  in  land  use  identification.  This  study  was  carried 
out  as  part  of  NASA’s  Water  Management  and  Control  ASVT  project, 

UCD  PROCEDURE 

An  operational  procedure  for  land  use  classification  from  LANDSAT  data 
has  been  developed  at  the  University  of  California,  Davis  (UCD)  for  use  by 
the  Corps  of  Engineers.  Referred  to  as  the  UCD  Procedure,  it  was  designed 
to  function  without  the  use  of  dedicated,  interactive  image  processing 
facilities.  Only  output  equipment  normally  available  in  Corps  field  offices 
(e.g.,  line  printer)  and  batch-mode  access  to  a  general  purpose  computer 
could  be  expected.  It  was  also  intended  that  the  procedure  would  not  re¬ 
quire  specialized  technical  expertise  in  data  analysis,  computer  program¬ 
ming,  or  remote  sensing. 
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The  UCD  Procedure  consists  of  an  organized  set  of  computer  programs  and 
manual  operations  for  the  identification  of  land  use  from  raw  LANDSAT  data. 

A  detailed  description  is  given  by  Algazi  (19791.  The  primary  tasks  of  the 
procedure  are: 

(1)  Obtain  LANDSAT  Computer-Compatible  Tapes  (CCT),  NASA  high  altitude 
ar‘rial  photography,  and  USGS  topographic  maps  for  the  location  and  date  of 
interest.  Extract  a  rectangular  area  of  data  containing  the  watershed  from 
the  CCT.  Check  for  radiometric  errors  in  the  LANDSAT  digital  data  and,  if 
necessary,  correct. 

(2)  Determine  the  geometric  registration  of  the  LANDSAT  image  with  the 
UTM  coordinate  system  of  the  USGS  topographic  maps,  LANDSAT  control  points 
are  identified  from  the  output  of  a  UCD  computer  program  which  enhances 
roads  and  water  bodies  found  in  the  LANDSAT  image.  A  regression  equation, 
calculated  from  the  two  sets  of  control  points,  provides  a  transformation 
mechanism  for  going  from  the  image  coordinate  system  to  the  map  coordinate 

s vs  tern. 


(3)  Use  an  unsupervised  clustering  algorithm  to  partition  the  LANDSAT 
foul — dimensional  data  space.  Groups  or  ’’clusters"  are  identified  that  con¬ 
tain  points  with  spectral  reflectance  values  that  are  simlar  to  members  of 
the  same  cluster,  and  dissimilar  to  the  points  of  other  clusters.  Each 
pixel  in  the  watershed  data  file  is  assigned  to  a  cluster. 
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(4)  Select  from  a  line  print(*r  map  of  the  cluster  assignments  six  sr*t  s 
of  adjacent  pixels  (spatial  groups all  belonging  to  the  samt'  ^'oster. 

Their  ct)rrespond  i  ng  location  on  th('  topograpliic  maps  is  determined  usin.g  the 
transformation  equation  of  step  (2),  Visual  translation,  from  the  map  to 
the  aerial  photographs,  of  the  spatial  group's  location  permits  a  land  use 
to  be  assigned  to  each  such  spatial  group.  For  clusters  having  a  consistent 
land  use  assigned  to  all  six  spatial  groups,  a  final  land  use  has  been  de¬ 
termined.  But  for  those  clusters  where  conflicts  exist  between  the  land  use 
identified  with  each  of  the  six  spatial  groups,  further  partitioning  of  the 
data  space  is  required. 

(5)  Clusters  with  conflicting  land  use  assignments  and  clusters  whose 
associated  land  use  could  not  be  determined  from  the  available  maps  and 
photos  are  reclustered  by  repeating  step  (3),  and  given  final  land  use  as¬ 
signments  by  repeating  step  (4). 

(6)  At  this  point  the  watershed  data  file  contains  a  land  use  classi¬ 
fication  (typically  5  to  7  categories)  for  all  pixels.  The  watershed  file 
is  then  resampled  at  the  grid  cell  centroids  using  a  nearest-neighbor  algo¬ 
rithm, 

(7)  The  resampled  file  is  then  entered  directly  into  the  watershed's 
grid  cell  data  bank  (as  explained  later).  Alternatively,  a  file  containing 
the  digitized  watershed  boundary  can  be  used  to  mask  the  resampled  file, 
leaving  only  the  grid  cells  within  the  boundary.  Total  acreage  of  each  land 
use  class  for  the  entire  watershed  is  then  computed. 
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I 'ANDSA'r -CON^E NTIONAL  LAND  USE  COMPARISON 

The  IJCD  procedure  lias  been  tested  at  HEC  on  two  watersheds:  Crow  Creek 
near  Davenport ,  Iowa  and  Walnut  Creek  near  Aus t in,  Texas ,  The  purpose  of 
the  applications  was  to  gain  familiarity  with  the  procedure  and  to  evaluate 
the  accuracy  of  LANDSAT-der i ved  land  use.  The  results  of  the  Crow  Creek 
classification  have  been  reported  by  Algazi  (1979). 

Walnut  Creek  has  a  drainage  area  of  55  square  miles.  LANDSAT  imagery 
for  3  May  1976  was  analyzed  using  January  1974  NASA  high  altitude  color  in¬ 
frared  aerial  photography  (scale  1:121,000)  and  photorevised  1973  USGS 
7-1/2-rainute  topographic  maps  (scale  1:24,000)  as  ground  truth.  Figure  1 
shows  the  location  of  the  watershed  on  the  LANDSAT  scene.  Six  land  use 
classes  were  identified:  cropland/pasture,  residential,  commercial/  indus¬ 
trial,  quarry,  forest/rangeland,  and  water. 

LANDSAT  land  use  was  entered  into  a  grid  cell  data  bank  that  had  prr- 
viously  been  constructed  by  the  Fort  Worth  District,  Corps  of  Engineers  for 
an  Expanded  Flood  Plain  Information  Study  (XFPI)  of  the  Walnut  Creek  basin. 
The  grid  cell  size  of  the  Walnut  Creek  data  bank  is  200  feet  (east-west)  by 
250  feet  (north-south),  or  1.148  acres.  Using  a  line  printer  spacing  of  8 
lines  per  inch,  line  printer  maps  of  this  grid  cell  size  correspond  to  the 
1:24,000  scale  of  USGS  7-1/2-minute  topographic  maps. 

Existing  land  use  had  been  determined  conventionally  from  manual  inter¬ 
pretation  of  October  1977  low  altitude  color  infrared  aerial  photography 
(scale  1:12,000).  This  land  use  classification,  referred  to  as  the 
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convent  i  on  a  \  land  use,  consisted  of  19  land  us(‘  cat(‘^ori**s  as  C-ivrn  in  'Iai>l«- 


1.  Table  1  also  lists  the  acrcMge  and  pe^rcent  distribution  for  t  bt' 

LANDSAT  land  use  categories.  Figure  2  contains  land  use  maps  for  botn  cia^^- 
s i f i cat i ons . 

A  cell -by-cel  I  comparison  of  the  ent  i  rt^  waterslied  required  that  the 
larger  number  of  conventional  land  use  categories  he  aggregated  to  the  f  t'V-  r 
LANDSAl  land  use  categories  as  shovsm  in  Table  1,  Because  conventional  ap(i 
lANBSAT  categories  are  not  always  compatible  this  created  some  problems.  As 
an  example,  the  conventional  category  ” t ranspor t a t i on/ commun i ca t i on/ 
utilities"  includes  major  highways,  right-of-ways  for  railroads  and  power 
transmission  lines,  airport  facilities  f including  buildings,  runwavs,  and 
vacant  land  within  the  airport  limits'),  and  sewage  treatment  facilities. 
LANDSAT,  however,  will  classify  the  water  surface  of  a  treatment  plant  as 
"water"  the  open  fields  surrounding  a  runwav  as  one  of  the  vegetation 
categories,  and  right-of-ways  as  whatever  land  surface  class  if  nearbv. 
Caution  should,  therefore,  be  used  in  interpreting  the  grid  cell  compar’son. 

Table  2  summarizes  the  cell-by-cell  comparison  of  LANDSAT  and  con¬ 
ventional  land  use  classifications.  Theoretically,  if  the  land  use  cate¬ 
gories  were  completely  consistent,  then  the  percent  of  grid  cells  that 
appear  down  the  diagonal  of  the  table,  52.7%,  would  represent  the  accuracy 
of  the  LANDSAT  c 1  ass i f i ca t i on .  An  additional  complication  is  the  inclusion 
of  "rangeland"  in  the  "cropland/pasture"  convent ional  land  use  category,  and 
in  the  "forest"  LANDSAT  land  use  category.  Being  unable  to  differentiate 
between  them,  one  should  consider  crop! and/pasture/forest/rangcl and  as  one 
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catogorv;  summing  the  diagonal  percentage's  would  then  ft'sul  t  in  70.8“1!  of  th 
grid  cells  "correctly*'  classified  by  LANDSAT . 

A  comparison  at  the  grid  cell  level  is  nearly  the  same'  as  a  comparison 
at  the  I^NDSAT  pixel  level;  both  Walnut  Creek  grid  cell  and  LANDSA7'  pixel 
are  approximately  1.1  acres.  At  svich  a  scale  errors  in  the  geometric  cor¬ 
rection  and  the  resampling  procedure  will  have  significant  impact  on  the 
computed  accuracy.  Another  comparison,  less  sensitive  to  such  factors,  is 
of  the  major  land  use  categories  at  tlie  watershed  level.  Looking  at  the 
total  percent  classified  as  (a)  residential,  Cb)  commerc i al / i ndus t r i al ,  and 
cropl and /pas tare/  f ore s t / range 1  and  by  conventional  and  LANDSAT  shows  a 
difference  of  7.6%,  and  -1.1%,  respectively;  the  average  absolute 

difference  for  the  major  land  use  categories  is  4.4%. 

LANDSAT  land  use  has  been  determined  by  UCD  staff  for  three  other 
basins:  Castro  Valley,  California;  Pennypack  Creek  in  Philadelphia, 

Pennysy 1 vania ;  and  Rowlett  Creek  near  Dallas,  Texas.  Corps  of  Engineers’ 
district  offices  had  determined  land  use  by  conventional  means  for  their  re 
spective  basins  and  encoded  this  information  into  grid  cell  data  banks. 
Cell-by-cell  comparisons  were  made  at  HEC  for  the  three  basins.  Results  of 
these  comparisons  were  consistent  with  the  findings  in  the  Walnut  Creek 
analysis:  individual  grid  cells  were  incorrectly  classified  by  LANDSAT 

approximate  1 y  38%  of  the  time,  whereas  aggregation  of  grid  cells  over  the 
entire  watershed  showed  mi  sc  1  a ssi f i ca t ion  of  the  major  land  use  categories 
averaging  2  to  8%, 


GRID  CELL  SPATIAL  DATA  MANAGEMENT  SYSTEM 


The  grid  cell  spatial  data  management  system,  with  which  the  foregoing 
maps  and  cell -by-cell  analyses  were  made,  has  been  tin  operational  tool  in 
the  Hydrologic  Engineering  Center  since  1975  ^U.S.  Army  Corps  of  Engineers, 
1975).  The  data  management  system  consists  of  a  set  of  utility  programs 
for:  fl)  encoding,  checking,  and  placing  the  geographic  map  data  in  the 

grid  cell  data  bank;  f2)  displaying  the  data  through  plotting  maps  of  one  or 
more  variables  by  their  absolute  values  or  by  weighted  combinations  of  rela¬ 
tive  attractiveness;  and  f3)  extracting  data  from  the  data  bank  and  formu¬ 
lating  parameters  for  hydrologic,  economic  and  environmental  simulation  pro¬ 
grams.  The  principal  computer  ptograms  for  the  LANDSAT  land  use  identi¬ 
fication,  data  bank  input,  hydrologic  parameter  identification  and  watershed 
model  are  shown  in  Figure  3. 

The  spatial  data  management  system  has  been  the  focal  point  of  the 
Corps  of  Engineers'  Expanded  Flood  Plain  Information  Studies  fDavis,  1978^. 
This  technology  was  developed  to  provide  a  comprehensive  management  tool  for 
use  by  local  governments  responsible  for  the  nation's  flood  plains.  The 
pilot  studies  undertaken  by  the  Corps'  district  offices  analyzed  the  hydro- 
logic,  economic,  and  environmental  aspects  of  existing  and  alternative 
future  land  use  patterns.  The  analyses  included  the  automatic  computation 
of  expected  annual  flood  damages  which  resulted  from  changed  hydrologic  re¬ 
sponses  and/or  stage-damage  functions  as  related  to  each  land  use  pattern 
and  flood  management  measure. 
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land  iis»'  i  dtMU  i  f  i  r  a  t  i  on  is  essential  to  the  foregoing 
analyses,  especially  for  the  economic  analysis  of  flood  damages.  Conse¬ 
quently,  twenty  or  more  separate  land  uses  were  often  required  to  accurately 
represent  the  flood  damage  relationships.  LANDSAT  was  found  to  be  capable 
of  identifying  about  six  land  uses  and  thus  could  not  be  used  for  the  eco¬ 
nomic  analyses.  The  LANDSAT  identified  land  use  was  thought  to  be  adequate 
for  hydrologic  purposes  and  the  ongoing  XFPI  studies  were  used  as  the  format 
to  compare  the  hydrologic  implications  of  lANDSAT  vs,  conventional  land  use 
classi ficati ons , 

The  HYDPAR  computer  program  (U.S.  Army  Corps  of  Engineers,  1978),  see 
Figure  3,  automatically  determines  hydrologic  parameters  from  the  geographic 
grid  cell  data  bank.  The  key  element  in  this  determination  is  an  input  re¬ 
lationship  between  the  geographic  features  and  the  hydrologic  parameters, 
e.g..  Table  3.  Any  rai nf a  1 1 /snowmel t- to-runo f f  simulation  methods  which  can 
be  related  to  geographic  features  can  be  handled  in  this  manner. 

HYDROLOGIC  COMPARISON 


The  primary  reason  for  examining  the  land  use  classification  ability  of 
LANDSAT  was  for  its  potential  application  to  hydrologic  modeling.  Cali¬ 
bration  of  hydrologic  models  typically  used  by  the  Corps  of  Engineers  in 
urban  areas  is  heavily  dependent  on  land  use  data,  particularly  in  basins 
where  land  use  is  changing  and  where  future  conditions  are  of  interest. 

The  computer  program  HEC-1  CU.S,  Army  Corps  of  Engineers,  1973)  has  the 
capability  of  explicitly  relating  land  use  to  runoff  using  two  procedures: 
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Snyder's  unit  hydrograph  with  percent  impt‘rvi  ousnt'ss ,  and  the  SCS  *  u:  v.‘ 
number  and  unit  hydrograph  ^U.S.  Soil  Conservation  Service,  1972'.  As  de¬ 
scribed  in  the  previous  section,  the  HYDPAR  program  accesses  the  nerrc^-jj-Y 
variables  from  a  grid  cell  data  bank  and  computes  the  specified  hydrologic 
parameters,  which  are  in  turn  input  into  an  HEC-1  model  of  the  basin. 

HYDPAR  contains  a  regression  equation  formulation  of  Snyder’s  lag  as  a 
function  of  stream  length,  length  to  centroid  of  subbasin,  stream  slop^,  and 
percent  imperviousness,  A  table  associating  a  percent  imperviousness  with 
each  land  use  category  in  the  data  bank  enables  HYDPAR  to  compute  subbasin 
percent  imperviousness  from  subbasin  land  use  distribution. 

In  a  similar  manner  HYDPAR  can  determine  the  SCS  unit  hydrograph  param¬ 
eter  from  stream  length,  basin  average  land  slope,  and  subbasin  average 
curve  number.  Curve  numbers  represent  an  empirical  relationship  betvjoen  hy¬ 
drologic  soil  type,  land  use,  and  their  resultant  runoff  potential.  From  a 
table  identifying  a  curve  number  with  each  combination  of  land  use  and  hy¬ 
drologic  soil  type,  Table  3,  HYDPAR  computes  subbasin  average  curve  number. 
Figure  A  illustrates  HYDPAR's  subbasin  printout  for  the  computation  of  the 
SCS  parameters. 

Both  procedures,  Snyder's  and  SCS,  wt re  used  te  l('st  the  "hvdrologic" 
accuracy  of  LANDSAT-der  i  ved  land  use.  Percent  i  mpt'rv  i  ousiiess  was  assigned 
to  each  of  Rowlett  Creek's  LANDSAT  and  conventional  land  use  categories. 
HYDPAR  was  used  to  compute  Snyder's  lag  for  a  portion  of  the  watershed 
f Upper  Spring  Cr.,  24. 6  square  miles)  and  a  cal i brat  i‘d  HEC- I  model  was  used 
to  simulate  runoff  from  selected  recurrence  interval  rainfall.  Differences 


11 


hotwot'n  the  discharge  frequency  curves,  based  on  either  I^NDSAT  or  conven- 
rional  land  use,  ran  be  interpreted  as  a  measure  of  the  hydrologic  signi- 
ficance  of  LANDSAT's  misclasslfication  of  land  use.  Considering  the  un* 
certainty  involved  in  estimating  any  frequency  curve,  the  difference  between 
Rowlett's  LANDSAT  and  conventional  curves  was  judged  insignificant.  In  the 
interest  of  conserving  space,  only  the  SCS  results,  described  in  the  next 
paragraph,  are  shown  graphically. 

The  SCS  method  was  applied  to  Pennypack  Creek.  The  SCS  curve  number 
and  lag  were  computed  by  HYDPAR  from  both  LANDSAT  and  conventional  land 
use.  Once  again  discharge  frequency  curves  were  estimated  from  the  sirnu- 
lated  runoff  of  a  calibrated  HEC-l  model,  this  time  using  SCS  parameters  to 
define  subbasin  rainfal l-runof f  response.  The  frequency  curves  of  Figure  5, 
based  on  LANDSAT  and  conventional  land  use,  are  nearly  identical.  Also 
shown  on  Figure  5  are  discharge  frequency  curves  which  would  result  if  the 
basin  were  modeled  by  all  industrial  and  all  natural  vegetation  curve  num¬ 
bers.  These  curves  demonstrate  the  possible  extremes  that  could  have  been 
generated  from  the  model  and  provide  a  set  of  references  from  which  to  eval¬ 
uate  the  actual  differences  between  LANDSAT  and  conventionally-based  fre¬ 
quency  curves. 
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CONCLUSIONS 

An  operational  procedure  for  determining  land  use  from  LANDSAT  imagorv 
has  been  applied  to  five  watersheds.  Based  on  our  experience  the  following 
conclusions  can  be  made: 

n)  At  the  grid  cell  level  LANDSAT  land  use  can  be  expected  to  be  in 
error  about  1/3  of  the  time. 

(2)  By  aggregating  the  land  use  over  the  entire  watershed,  LANDSAT' s 
misclassi fication  of  land  use  reduces  to  2  to  8%  for  the  major  land  use  cat- 
egor i es. 

C3)  Both  of  the  above  accuracy  assessments  must  be  qualified  by  noting 
that  conventional  land  use,  to  which  LANDSAT  was  compared,  will  sometimes 
have  land  use  categories  that  are  inconsistent  with  the  LANDSAT  land  use 
categories  (i.e.,  we  are  sometimes  comparing  apples  and  oranges'^.  Also, 
errors  introduced  during  geometric  correction  and  resampling  will  he  inter¬ 
preted  as  LANDSAT  misclassification  errors  when  comparing  on  a  cell -by-cell 
basis. 

(4)  The  UCD  Procedure  works.  It  is  a  complete,  self-contained  package 
of  computer  programs  and  manual  operations  that  permit  a  user  to  identify 
land  use  from  lANDSAT  digital  data  without  requiring  the  use  of  expensive, 
interactive  image  processing  equipment. 
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f  S  )  iw.'i  1  iiat  tnl  in  tf^rms  of  tho  riifforrnce  in  di  scharge  froqnoncy  curvf^s, 
the  LANDSAT-der  i  ved  land  use  was  found  to  be  completely  adequate*,  Th.*  num¬ 
ber  and  type  of  land  use  categories  derived  from  LANDSAT  data  W('re  suf¬ 
ficient  to  be  able  to  apply  two  standard  ftydr^^logic  modeling  techniques, 
Snyder's  unit  hydrograph  with  percent  imporvi ousness  and  the  SCS  curve  n tim¬ 
ber  and  lag  method, 

(b)  LANDS AT  land  use  can  be  directly  incorporated  into  a  watershed's 
grid  cell  data  bank,  thus  providing  an  automated  environment  for  applying 
the  IJVNDSAT  classification  in  routine  hydrologic  investigations. 
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tabu;  1 

WAIJ4UT  CREEK  LAND  USE  DISTRIBUTION 
CONVENTIONAL  AND  lANDSAT 


Convent i ona 1  Land  Use 

Acres 

Z 

LANDSAT  Land  Use 

Acres 

c. 

Low  density  residential 

1,785 

4.9 

Resi dent i a  1 

9,135 

25.0 

Med  density  residential 

3,671 

10.0 

Hi  density  residential 

156 

0.4 

Multi-farm  resi.  3ntial 

619 

1.7 

Mobile  homes 

126 

0.3 

Strip  commercial 

264 

0.7 

Commerc ial/industrial 

1,861 

5.1 

Shopping  centers 

85 

0.2 

Inst i tut i  na 1 

536 

1.5 

Indust rj  .i. 

488 

1.3 

Ind/com  ct  Tiplexes 

742 

2.0 

Pub  lie  iij 

96 

0.3 

Trans  f  comun/u t  i  1  i  t  ies 

1,261 

3.4 

Barren  land/quarry 

1,022 

2.8 

1 

Barren  land /quarry 

i 

223 

_ ! 

0.6 

Cropland 

i 

3,917 

10.7 

Cropland/ pasture 

15,298  ! 

41.8 

Pasture/ rangeland 

11,327 

31.0 

Dev  open  space 

139 

0.4 

Undev  open  space 

1,135 

'  3.1 

Forest 

'  9,143 

"  25.0 

Forest/ rangeland 

9,988 

27.3 

Water 

62 

0.2 

Water 

70 

! 

0.2 

36,574 

99.9 

36,575 

TABLE  3 


PENNYPACK  CREEK  LANDSAT  LAND  USE 

CURVE  NUMBER  SUMMARY  WITH 
ASSOCIATED  LAND  USE  CATEGORIES 


LAND  USE 

HYDROLOGIC 

SOIL 

TYPE 

CATEGORY 

TITLE 

A 

B 

c 

D 

1 

Natural  Vegetation 

40. 

63. 

75. 

81. 

2 

Pasture 

62. 

75. 

83. 

84. 

3 

Agricul tural 

62. 

75, 

83, 

84. 

4 

Developed  Open  Space 

39. 

61. 

74. 

80. 

5 

Residential 

65. 

78. 

85. 

88. 

6 

Light  Industry 

81. 

88. 

91. 

93. 

7 

Intermediate  Industry 

81. 

88. 

91. 

93. 

8 

Heavy  Industry 

81. 

88. 

91. 

93. 

9 

Water 

100. 

100. 

100. 

100. 

Figure  2. 

WALNUT  CREEK  LAND  USE  COMPARISON  MAPS 
(Darker  to  lighter  shades  represent  Industrial/Commercial ,  Residential, 
Natural  Vegetation,  and  Agricultural  land  use  categories,  respectively.) 
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Figure  3.  SCHEMATIC  OF  PROCESSING  SYSTEM 


tXCrtOANCE  FREQUENCY,  in  percent 


Figure  5 

PENNYPACK  CREEK  BASIN 
SUBBASINS  1-65  (55.74  sq.mL) 
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